PROCEE DINGS 


AMERICAN SOCIETY 


ENGINEERS 


JULY, 1952 


DISCUSSION 
SURFACE CURVES FOR STEADY 
NONUNIFORM FLOW 
(Published October, 1951) 


Ivan Nelidov; Stevens; Maurice Bishay 
and Robert Jansen 


HYDRAULICS DIVISION 


Copyright 1952 the AMERICAN ENGINEERS 
Printed the United States America 


Headquarters the Society 
St. 
New York 18, N.Y. 


PRICE $0.50 PER COPY 


SOCIETY OF 
CIVIL 
ENGINEERS 


GUIDEPOST FOR TECHNICAL READERS 


value significance readers various fields are here listed, 
for convenience, terms the Society’s Technical Divisions. Where there seems 
overlapping interest between Divisions, the same Separate number may appear under 
more than one item. 


Technical Division Proceedings-Separate Number 
42, 43, 48, 52, 60, 93, 94, 95, 100, 103, 104, 108, 121, 130 
(Discussion: D-X XVIII, D-7, D-16, D-18, D-23, D-43, 
D-75) 
58, 60, 62, 64, 93, 94. 99, 101, 104, 105, 115, 131, 138 (Dis- 
cussion: D-16, D-23, D-43, D-60, D-62, D-65, D-86) 
43, 50, 55, 71, 92, 94, 103, 108, 109, 113, 117, 121, 126, 


130, 132, 133, 136, 137 (Discussion: D-3, D-8, D-17, 
D-23, D-36, D-75, D-92) 


Engineering Mechanics............ 122, 124, 125, 126, 127, 128, 129, 134, 135, 136, 139, 141 
(Discussion: D-24, D-33, D-34, D-49, D-54, D-61, D-96, 
D-100) 

43, 44, 48, 58, 70, 100, 105, 108, 113, 120, 121, 130, 137, 
138 (Discussion: D-X XVIII, D-23, D-60, D-75) 

50, 55, 56, 57, 70, 71, 78, 79, 80, 83, 86, 92, 96, 106, 


107, 110, 111, 112, 113, 116, 120, 123, 130, 134, 135, 139, 
141 (Discussion: D-70, D-71, D-76, D-78, D-79, D-86, 
D-92, D-96) 

Irrigation and Drainage........... 97, 98, 99, 102, 106, 109, 110, 111, 112, 114, 117, 118, 120, 
129, 130, 133, 134, 135, 138, 139, 140, 141 (Discussion: 
D-XXIII, D-3, D-7, D-11, D-17, D-19, D-25-K, D-29, 
D-30, D-38, D-40, D-44, D-47, D-57, D-70, D-71, D-76, 
D-78, D-80, D-86, D-87, D-92, D-96) 

48, 55, 56, 69, 71, 88, 96, 103, 106, 109, 110, 117, 118, 
120, 129, 130, 133, 134, 135, 139, 141 (Discussion: 
D-XXIII, D-2, D-3, D-7, D-38, D-40, D-44, D-70, D-71, 
D-76, D-78, D-79, D-86, D-92, D-96) 

Sanitary Engineering............. 55, 56, 87, 91, 96, 106, 111, 118, 130, 133, 134, 135, 139, 
141 (Discussion: D-29, D-37, D-56, D-60, D-70, D-76, 
D-79, D-80, D-84, D-86, D-87, D-92, D-96) 

Soil Mechanics and 44, 48, 94, 102, 103, 106, 108, 109, 115, 130 (Discus- 
sion: D-4, D-XXVIII, D-7, D-43, D-44, D-56, D-75 
D-86) 

42, 49, 51, 53, 54, 59, 61, 66, 89, 100, 103, 109, 113, 116, 
117, 119, 121, 122, 123, 124, 125, 126, 127, 128, 129, 132, 
133, 136, 137 (Discussion: D-51, D-53, D-54, D-59, D-61, 
D-66, D-72, D-100) 

Surveying and Mapping........... 50, 52, 55, 60, 63, 65, 68, 121, 138 (Discussion: D-60, D-65) 

41, 44, 45, 50, 56, 57, 70, 71, 96, 107, 112, 113, 115, 
120, 123, 130, 135 (Discussion:, D-8, D-9, D-19, D-27, 
D-28, D-56, D-70, D-71, D-78, D-79, D-80) 


constant effort made supply technical material Society members, over the entire 
range possible interest. Insofar your specialty may covered inadequately the 
foregoing list, this fact gage the need for your help toward improvement. Those who 
are planning papers for submission will expedite Division and 
Committee action measurably first studying the ASCE for Development 
Proceedings-Separates” style, content, and format. For copy this Manual, address 
the Manager, Technical Publications, ASCE, 39th Street, New York 18, 


The not responsible for any statement made opinion expressed 
tts publications 
Published Prince and Lemon Streets, Lancaster, Pa., the American Society 
Civil Engineers. Editorial and General Offices West Thirty-ninth Street, 
New York 18, Reprints from this publication may made 
condition that the full title paper, name author, page 
reference, and date publication the Society are given. 


Proceedings—Vol. July, 1952 Separate No. D-96 
DISCUSSION 


Ivan sketching the water-surface curves obtained with 
the aid the slopes water, computed means Eq. 12, the author has 
made the water-slope tangents intersect the midpoints the depth incre- 
ments. This equation may rewritten follows: 


which (conforming the notation the paper), the friction slope and 
the specific energy, energy referred the invert. Eq. can repre- 
sented the Bernouilli equation the depth added both sides, the 
friction slope introduced, and both sides are multiplied dz; thus, 


which the subscripts and refer two sections, distance apart; 
the difference elevation between the inverts section and 
section and the friction head lost. Eq. may also re- 


written as, 


which and are the specific energies sections and respectively. 

The only correct way construct the water-surface curve integration. 
However, cases where integration can used are rather limited because the 
channel not prism uniformly varying cross sections. may quite 
irregular and can change abruptly from section section. 

The approximation proposed the author involves constructing the water- 
surface profile means tangents conneeting two consecutive depths such 
way that these tangents meet the midpoint the depth increment. 
may noted, that the standard mathematical procedure for constructing 
curve mark equal horizontal increments dz, with varying vertical incre- 
ments dy. 

The tangents will thus intersect the vertical ordinates or, the ordinates 
are translated distance they will intersect the midpoint between the 
ordinates. 


paper Robert Jansen was published October, 1951, Proceedings-Separate 
No. 96. The numbering footnotes, equations, illustrations, and tables this Separate are continua- 
tion the consecutive numbering used the original paper. 


Civ. and Structural Engr., San Francisco, Calif. 
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example the ordinates 2.5 and 3.0 are selected, the horizontal dis- 


0.25 0.25 

0.00275 0.00303 173.4 ft. The same 
computed the Bakhmeteff method,? and using his notation, obtained 
follows: 6.12 ft; 0.001; the hydraulic exponent 2.98; (the 


critical slope) 0.00244; 0.41; and 


tance between them will be: 


6.12 


which the varied flow functions are 0.5059 and 0.4155 
The curve this type and the flow rapid and retarded. 
use the Bernoulli equation determining distance without trial, Eq. 
must solved, from which: 


The numerical values entering Eq. are follows: 2.50 7.25 9.75; 
basis average hydraulic radius Ray 2.155 and average velocity 


9:75 8.02 1.73 
19.8 per sec, with 0.013. Then, 0.01095 0.001 0.00995 
173.6 ft. 


average value the friction slopes computed for the end section used, 
the distance would differ more from the foregoing values. The following 
numerical data would used: 21.6 per sec; 2.00 ft; 18.0 
per sec; and 2.31 ft. Furthermore, 0.01455, 0.00815, and 
loss should computed introducing Ray and Vay for end sections and 
determining the friction slope the basis their values. This may proved 
taking the general expression for friction loss and comparing the results. 
Thus, integrating the theoretical friction loss with the aid the Bakhmeteff 


yi)? 
1932. 
76. 
82. 


a 2 
dit = (16) | 
| 
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being the conveyance the normal depth and, the distance along 
the line flow from initial section. 

friction-loss computation based the data referred the average 
cross section, the expression for friction loss will be: 


3 


loss computed taking the average friction slope, the 
expression then will 


3 


Table indicates that friction computations the Bernoulli equation 
should based average area and average hydraulic radius, and not 
average friction slope. general, the method proposed the author gives 
results within the desired degree accuracy—that is, within the ordinary 


VALUES 
Friction 


Equation 
0.7 1.3 
Constant [3.42] Constant [1.36 
Constant [3.93], +15% Constant [1.45], +7% 


range variation depths. However, with rapidly varying depths, all 
approximation methods (which includes the author’s method and the Bernoulli 
method when finite distances are used) are likely give incorrect results, 
the degree error being established only after integration performed. 

important classify the water-surface curves properly that the 
type involved evident glance. Fig. three major classes are noted, 
first: flow with (class A), with (class B), and with 
(class C). 

The slope invert class adverse (class zero (class and 
gentle (class class critical; and class greater than critical. 
any horizontal row the first figure shows the water depth greater than both 


» 
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(m) TRANQUIL, RETARDED RAPID, ACCELERATED RAPID, RETARDED 
C.—STEEP SLOPE, 


WATER SURFACE CuRVES STEADY NONUNIFORM 


None 
f (a) (6) TRANQUIL, ACCELERATED (c) RAPID, RETARDED 
SLOPE, 
ay <0 dz ‘ 
(d) (e) TRANQUIL, ACCELERATED (f) RAPID, RETARDED | 
BOTTOM, 
>0 a <0 
(g) TRANQUIL, RETARDED TRANQUIL, ACCELERATED RAPID, RETARDED 
SLOPE, 
y 
dz 
Wg 
TRANQUIL, RETARDED (k) UNIFORM RAPID, RETARDED 
SLOPE, 
dE 
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and the second figure between and and the third figure 
smaller than both and 

Consistent with the letter symbols adopted for the author, the notation for 
water-surface curves the same that introduced Mr. Bakhmeteff and 
the classification itself only more graphical presentation Mr. Bakhmeteff’s 
ideas. far the type the flow concerned, the first vertical column 
the flow tranquil and retarded; the second vertical column tranquil 
and accelerated until becomes flow critical slope, where uniform, and 
past it, becomes rapid and accelerated; the third vertical column the flow, 
for all cases rapid and retarded. 

The writer has found that this form classification helpful determining 
the type curve involved before beginning the detailed computations. 


really simple problem the hardest possible way. Since 1925, the writer has 
been urging the use simple step-by-step method applicable all 
This method not original with the writer. was first presented Alva 
strictly applicable any form water surface profile 
sustaining adverse slopes under sub-critical super-critical flow rough 
smooth channels provided only they are 
uniform channels—that is, rectangular, tri- 
angular, trapezoidal, circular, parabolic, ellipti- 
cal, otherwise—just that the area and, 
consequently, the perimeter, radii, are the 
same functions the depth through the reach 
under 

illustrate the point that the writer 
wishes make will necessary compute 
the surface curves Example the Husted 
method. 

Fig. illustrates the conditions for this ex- 
ample which the bed slope greater than 
the friction slope, the depth super-critical and super-normal, and there 
gain specific energy downstream. From Fig. 


Friction Slope 


and 
(20c) 
which 


Koon, Cons. Engrs., National Bldg., Portland, Ore. 
Engineering News-Record, Vol. 95, 1925, 550. 
Transactions, ASCE, Vol. 102, 1937, 666. 


Method Computing Backwater and Drop-Down Curves,” Alva Husted, Engineering 
News-Record, Vol. 92, 1924, 719, 
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the specific energy water measured head above stream bed 
datum. Dimension finite length channel over which the surface 
the stream considered straight line; the bed slope; and the 
friction slope. Table shows the computations for this example. The flow 


TABLE FoR EXAMPLE THE METHOD 


slope; and Cols. and 13, finite length channel, and Az, the distance from the starting 
point the surface curve. 


(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 


Hydraulics,” Horace King, McGraw-Hill Book Co., Inc., New York, Y., 
Ed., 1939, 309. 


the direction increasing depths. There gain specific energy down- 
stream shown Fig. due the accelerating force shows 
the surface curve taken from Fig. circles are points given Table 
The agreement practically perfect. 


800 1000 1200 1400 1800 2000 
Distance Along Channel Feet 


Depth Water (y) Feet 


Fia. 


Example the same that presented Mr. Bakhmeteff? 1932, who 
recognized that this surface curve includes hydraulic jump but did not define 
it. 

For this example the bed slope less than the friction slope, the depths are 
sub-critical and sub-normal, and specific energy diminishes downstream, result- 
ing the formula, 


4 
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Fig. shows the actual surface curve. The computations began the 
critical depth 4.49 and extended upstream using depth decrements 0.10 
account the rapidly changing curvature. The circles are points given 
Mr. Bakhmeteff. 

locate the position the hydraulic jump was necessary compute 
the surface curve downstream from the critical depth where super-critical 
depths obtain. points upstream from the critical depth the values 


were computed from both the sub-critical and super-critical curves, 


where the pressure the channel cross section and the 


momentum the water, set forth Julian ASCE. 


800 600 400 200 200 400 600 1000 1200 1400 1600 1800 2000 
Distance Along Channel, Feet 


The height the jump was determined the method given the 
4.24 above channel bed. The super-critical surface curve was then 
moved bodily upstream from point point The solid lines Fig. 
therefore represent the actual surface curve for the author’s Example neglect- 
ing energy losses. 

The computation surface curves from varied-flow function, used 
the author, Messrs. Bresse, Bakhmeteff, and others, practically precludes the 
possibility taking all the factors into account. the other hand, they are 
readily applied the step-by-step method. For example, the writer were 
trying find the surface curve for Example would compute from 
equation something like 


which equals the velocity head the mean velocity. 


Hydraulic Jump and Critical Depths the Design Hydraulic Julian Hinds, 
Engineering News-Record, Vol. 85, 1920, 1034. 
‘The Hydraulic Jump Standard Conduits, Stevens, Civil Engineering, October, 1933, 


700 
: g - 
Flow 
wu -— — 
i=) 3 | 
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much closer approximation the true surface curve Example from 
the sluice gate the foot the hydraulic jump would thus obtained, that is, 
adding the kinetic energy the mean velocity and then assuming 
only 80% velocity head recovery expanding areas. 


lem involving the computation backwater curve has been solved the 
author’s method and five other methods. Results this computation are 
shown Table and Fig. Comments are given the relative accuracy 
the various methods and the relative tediousness the work involved 
each method. 


7.0 
Solution 
Step 
6.5 
Bakhmeteff 
Ming Lee 


Distance Channel, Thousand Feet 


Elevation, Feet Above Datum 


Problem.—The problem selected for solution requires the construction the 
backwater curve for stream flowing under sluice gate into rectangular 
channel whose bottom width ft. The slope the bottom the channel 
0.0004 and the value Manning’s 0.01486. gate regulated 
that the rate discharge 136 per sec, with depth the vena contracta 
7.0 ft. 


Graduate Student, Univ. Illinois, Urbana, ill. 
Prof. San Eng., Univ. Illinois, Urbana, 
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Preliminary Steps depth uniform flow, computed 
the Manning formula, found Critical depth found 1.79 ft. 
Hence, the surface curve will the M-1 type gradually varied flow 


rectangular channel. 


explained ALONG THE CHANNEL, 
ASCE, but using the Ameri- 

Depth Graph- Bakh- Mono- 


hydraulic symbols (ASA 
5.0 
5.5 


Percentage deviation from the graphcial method depth 


which and represent any two points the backwater curve. That is, 
the distance between two points the backwater curve can obtained the 


8 


integration between points and The value can ob- 
tf 


tained from the expression, 


TABLE THE BACKWATER the curve plotted 


GRAPHICAL INTEGRATION against any con- 


venient scale, the distance 


(ft) (sq ft) ackwater curve 


(1) (2) (3) (4) (5) (6) computed 
4.5 2.37 4,643 2.92 9.26 
5.0 2.50 2.18 5.50 3,281 3,281 the area under the 
5.5 2.62 1.70 4.35 5,683 
6.0 2.72 6,080 1.35 3.775 1,955 7,638 
6.5 2.82 6,568 1.10 3.45 1,748 9,386 
7.0 2.92 7,059 0.902 3.23 1,639 11,025 -curve. Such 


curve has been drawn and 
the results the compu- 
tation graphical 
method areshownin Table6. 

Solution Step expression the distance between two points 
backwater curve can derived from the Bernouilli theorem and expressed 


along the channel below depth 4.5 ft. 


Mechanics Hunter Rouse, John Wiley Sons, Inc., New York, Y., 
1946, pp. 80-82. 
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standard nomenclature 


Results the computation the backwater curve this expression applied 
the illustrative problem are 


gration, Using the Bakhmeteff 


Specific 


4.5 4,643 2.92 


which can solved the 
computed Mr. Bakh- 

meteff, proceeding small 
increments depth that may assumed constant within each inter- 
val. Results computations Mr. Bakhmeteff’s method are compiled 
Tables and 

The Mononobe this method the assumption made that both 
the area and the wetted perimeter may expressed monomial func- 
tions the depth For example, equal constant times and 
equal constant times 
Mononobe, based the RECTANGULAR CHANNEL 


preceding assumptions, 
Wetted Hydrau- 


(1) (2) (3) (4) (6) 
which 
1.001 1.001 
and 
(31) 


Open Bakhmeteff, McGraw-Hill Book Co., Inc., New York, 
Y., 1932, Chapter VIII. 

Investigation the Backwater Profile for Steady Flow Prismatic Lans- 
ford and Mitchell, Bulletin Series No. 381, Univ. Illinois, Urbana, 1949. 

and Drop-Down Curves for Uniform Nagaho Mononobe, Transactions, 
ASCE, Vol. 103, 1938, 950. 


Ae, 
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The work involved the computation the complete backwater curve for 
the illustrative example tedious that only the distance between the depths 
4.5 and 7.0 has been located. was found means the equation 


9.—COMPUTATION THE BACKWATER CURVE 
THE BAKHMETEFF METHOD 


(1) (2) (3) (4) (6) (7) (8) 
4.5 1.125 0.125 0.9121 0.8045 
5.0 1.250 0.125 0.9146 0.9133 0.5760 —0.2285 3,336 3,336 
5.5 1.375 0.125 0.9174 0.9160 5,723 
6.0 1.500 0.125 0.9204 0.9190 0.3720 7,708 
6.5 1.625 0.125 0.9229 0.9216 0.3165 9,469 
7.0 1.750 0.125 0.9254 0.9242 0.264 11,204 
Distances feet, along the channel, below depth 4.5 ft. 
and charts that 
4.0 (3.4)? (2) 
Method Ming Lee’s basic expression the form, 
y2 N A? 


(1) (2) (3) (4) (5) (6) (7) 
4.5 0.790 1.308 
5.0 0.125 0.640 0.884 3,344 3,344 
5.5 0.125 0.530 0.667 5,737 
6.0 0.125 0.445 0.527 1,987 
6.5 0.125 —0.0555 0.378 0.427 9,483 
7.0 0.125 —0.0525 0.330 0.363 11,230 


Distances feet along the channel below depth =4.5 ft. 


which are given Mr. Using Mr. Bakhmeteff’s and Mr. Lee’s 
tables, values the coordinates points the curve have been computed, 
with the results shown Table 10. 


“Steady, Gradually Varied Flow Uniform Channels Mild Slopes” Ming Lee, thesis presented 
the University Illinois, Urbana, 1947, partial fulfilment the requirements for the degree 
Doctor Philosophy. 


The Jansen computations the slopes the curve 
curves this method the tangents are made intersect the midpoint 
between contiguous depths. The backwater curve resulting shown Fig. 

Summary.—Theoretically, the graphical method yields the most accurate 
results because assumptions are made concerning the conditions the 
integration procedure, other than those made the derivation the equation 
gradually varied flow. However, the work involved this method tedious. 
The work involved the step method relatively simple and the results are 
sufficient accuracy for most purposes. The methods direct integration pro- 
posed Mr. Bakhmeteff and Mr. Lee give similar results and are less 
tedious than the graphical method, but may less accurate because assump- 
tions made the computation the integral functions. the Bakhmeteff 

method, the empirical relationship 


BACKWATER CURVE THE depth holds for wide range 
JANSEN METHOD depths. However, there may 
appreciable changes the velocity 

Depth, Factor, Factor, Slope 
heads that are significant. Hence, 


the limitations imposed assuming 


(1) (2) (3) 


4.5 15.60 11.0 
5.0 11.82 0.000177 restrict the computations small 

7.29 4.635 0.000260 increments depth. Mr. Lee’s 


the depth, for practical purposes, 

that the integration limits may 
taken over wider ranges depth. Under such conditions the Lee method may 
convenient use. 

the Mononobe method the computations are more direct that the 
distance along the channel between any two sections may determined one 
step. The integration the gradually varied flow equation simplified 
making the assumptions that both the wetted perimeter and the area can 
expressed monomial functions the depth. The drawback the Mono- 
nobe method the difficulty using the graphical charts prepared him, and 
the inherent errors resulting from interpolations the reading the chart. 

Mr. Jansen’s method and Mr. Mononobe’s method differ from each other 
and the other four methods investigated, the results obtained the illustra- 
tive problem. The difficulties with these two methods lie mainly the reading 
the intersection lines and curves. the Jansen method, order ob- 
tain uniform curve the computed tangents are made intersect the mid- 
point between contiguous depths. The accuracy with which the location 
such points determined fixes the degree approximation the shape the 
resulting curve. 


open channels carrying steady nonuniform flow was first solved 


Lecturer Civ. Eng., Univ. Hong Kong, Hong Kong. 


for infinitely wide rectangular channel early 1860. 1912 
Mr. Bakhmeteff first suggested generalized solution for channels any shape 
which, 1932, published? giving many numerical solutions. 1948 
Phillip ASCE, proposed dimensionless formula with set 
curves for the solution the hydraulic exponent, ng, for rectangular and 
trapezoidal channels. 1949 Von ASCE, introduced 
another exponent order facilitate the integration the equation 
steady nonuniform flow. Mr. Jansen’s paper further extension, and, the 
same time, new and ingenious approach the same problem. stated 
the author, the method applicable only rectangular and trapezoidal 
channels symmetrical shapes. 

10, 11, and can extended that values the u-functions and 


j-functions are dimensionless forms. Substituting y/b and simplifying, 


Eq. becomes 


and becomes 


Next, assume that andj (both which are dimensionless 
then become 


and 
1+2¢a 
Eq. can then written 
which 
2 
(36a) 


Sons, Inc., New York, 

Kirpich, Civil Engineering, October, 1948, 47. 

the Equations Nonuniform Von Seggern, Transactions, ASCE, Vol. 
115, 1950, 71. 
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and 


All terms the extreme right side Eq. are dimensionless. For any 
prismatic trapezoidal channel (including rectangular channels), both and 


Values 


0.03 


0.01 
0.01 0.1 10? 104 105 
Values Functions 


can easily calculated, and and can obtained either from curves (such 
Fig. from prepared tables; and the slope between the desired 
sections can computed and 36. 


equation for the slope the water surface 
trapezoidal channel for steady nonuniform flow has been developed the 
author. The practical solution the equation requires the use special 
tables. The surface curve itself must found from various values the slope 
either graphically numerically. order test the advantage the pro- 
posed method, seems logical compare with the usual tabulation method. 


Director Water Section, Public Works Dept., Tel-Aviv, State Israel. 


ts 
/ 
2.0 
| 
| 
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The basic relation (Fig. 10) expressed 


With and known, textbooks usually assume definite value for 
find values and that satisfy Eq. 37. The method thus becomes subject 
trial-and-error, which can avoided solving Eq. for dz. inspection 
the right side Eq. shows that given, 
shooting flow exists. determine as- 
channel bottom, given and the energy 
gradient can computed the relation 


assuming that 


72 2 


The application the method shown Table 12, using Example 

The values are converted approximately into metric system, that is, 
culation was made slide rule. 

The result checks sufficiently with Fig. The tabulation method can 
used for any uniform cross section and gives direct solution for the water sur- 
face whereas the proposal the author restricted special shape cross 
section, requires previous tabulations, and includes (for the graphical solution) 
flat and, therefore, inaccurate intersections lines. The tabulation method, 
therefore, seems simpler and most easily visualized. 


TABLE SURFACE SLOPE, EXAMPLE 


(1) (2) (3) (7) (8) (9) (10) (11) (12) (13) 


1.1 6.75 0.60 8.35 4.15 2.60 96.3 0.44 0.33 0.0495 0.0485 53.5 
0.60 4.15 0.75 6.65 3.00 1.15 0.54 0.44 0.0220 0.0210 55.0 
6.65 3.00 0.90 5.55 2.47 0.53 37.5 0.64 0.55 0.0115 0.0105 50.5 


min 2.23693 miles per hr; and meter raised the four-thirds power (Col. 10) 133.8 


the case the more complicated problems actually encountered prac- 
tice, the idea finding the surface line calculating the slope first, can used 
with advantage tabulation method. For example, the writer has been re- 
quired calculate the backwater profile for trapezoidal river channel with 


/ 
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gradually increasing width and dikes for double profile. The water-filled 
cross section that case not function the depth alone. also neces- 
sary know the place along the channel which this depth dis- 
tance between control sections assumed and, for each these sections, the 
elevation and slope the energy grade line determined, which will result 
continuous line without sters the intersections. 


1:1250;Q per sec (263 acre-ft per 0.025; the flow conditions 
section are assumed known, with 200 and 1:1,000. The 
discharge conditions control sections are examined over distance 
200 with the energy grade lines interescting thus: 


Symbol Computation Metric units English units 
2.0 (22.0 4.0) square meters 560 square feet 
0.025 
1600 
1.68 meters inches 
1/1000 1/1000 
1/31.6 1/31.6 
1.78 meters 351 feet per 
31.6 
per second minute 
1.78 cubic meters 272 acre-feet 
per second per hour 


lon2 
Oo zx 


Symbol 
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Computation 


0.16 2.00 


100 
1000 


2.16 0.10 


100 
1250 


2.26 0.08 


Metric units 


0.16 meter 


2.00 meters 
2.16 meters 


0.10 meter 
2.26 meters 
0.08 meter 


2.18 meters 
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English units 


6.3 inches 


78.7 inches 
85.0 inches 


3.9 inches 
88.9 inches 
3.1 inches 


85.8 inches 


assumed that 2.03 this case, the correct value found 
most easily trial and error; but possible calculate, beforehand, 
table diagram several values from which that one selected which will 


fulfill the required conditions. 
thus (see Fig. 12): 


The influence the side channels neglected, 


Symbols Computations Metric units English units 
2.03 (19.4 4.0) 47.6 square meters 512 square feet 
0.03 4.0 
1.67 meters 65.7 inches 
1.99 meters 78.4 inches 
1.95 meters 384 feet per minute 
47.6 
per second 
3.82 
0.19 meter 7.5 inches 
2.03 80.0 inches 
0.19 2.03 2.22 meters 87.5 inches 
830 


| 
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Symbols Computations Metric units English units 
0.16 meter 6.3 inches 
1250 
2.22 0.16 2.38 meters 93.8 inches 
100 

330 0.12 meter 4.7 inches 
2.38 0.12 2.26 meters 89.1 inches 
2.38 0.12 2.50 meters 98.4 inches 
2.50 0.24 2.26 meters 89.1 inches 


The value 2.26 meters (89 in.) agrees with that found the preced- 
ing computation. 

With the continuation the calculation sections upstream will 
necessary take the flow the side channels into account. will suffi- 
cient distribute the discharge between the main channel and the side chan- 
nels accordance with the relation, 


and then continue the calculation for the main channel only. 

Separation losses and eddy losses can accounted for suitable selection 
the friction value for each stretch dz, separately. 

summarize: The basic idea offered this paper should not wasted 
simple problems for which the existing methods are good enough, but should 
used for more complicated cases, where the existing trial-and-error methods 
are inadequate. 


Rosert J.M. ASCE.—A lucid verification the results 
yielded the suggested graphical procedure has been presented Mr. 
Stevens. should not surprising that curves plotted this method con- 
form with those computed the equation Mr. Husted, since both methods 
solution are based the classic law the conservation energy, and both 
make use weighted energy gradients define the water slope finite in- 
crement channel. Any discrepany observed between the two curves should 
primarily reflection differences assumptions energy loss and 
channel increment. 


takes the form: 


25Civ. Engr. Associate, Bureau Eng., City Los Angeles, Calif. 
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Eq. can transformed 
(43) 


which the same Eqs. and 22. 

deriving Eq. 43, friction was the only cause energy loss included. 
order take eddy loss also into account, the term should sub- 
tracted from the right side Eq. that the final equation becomes 


lamentable that relationships for energy loss open channel have 
not been developed for nonuniform flow. The use the Chezy equation 
this case pardonable only because satisfactory substitute 
though generally accepted that the loss diverging flow greater than 
that uniform flow, reliable mathematical statement the difference 
available the technical literature. The application the eddy loss 
cient, isa rough attempt compensate for the lack Mr. 
Stevens’ selection, 0.20, value frequently applied practice. 

inject eddy-loss factor such this into the writer’s equations re- 
quires only simple modification. the velocity head differential Eq. 
multiplied the quantity suggested, can shown that Eq. 
will become 


Q? 
(45) 
and Eq. will changed 


Since both the coefficient and the factor applied Mr. Stevens, are 
dimensionless, their introduction into the writer’s equations calls for very 
minimum effort. This contradictory Mr. Stevens’ opinion that 
computation surface curves from varied-flow function used the 
author practically precludes the possibility taking all the factors into 

The writer does not consider necessary here defend the merits the 
Manning formula, particularly since, this case, yields results close those 
given the Kutter relation. However, Mr. Stevens may find some argument 
from other quarters his assertion that the Manning expression dimension- 
ally incorrect since the dimension the roughness factor topic for debate. 
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Contrary the opinion Mr. Stevens and Mr. Pietrkowski, the tables 
the depth functions, and are not essential the method presented the 
paper—no more so, fact, than they would the Husted procedure 
methods advocated others. will noted that these functions are simply 


and 
(48) 


The Manning value the energy gradient obtained multiplying 


the term thus, 


The significance the function defining the change velocity head has 
already been indicated Eq. That relation can rewritten give 


apparent from Eqs. and that the application the depth func- 
tions and the computation surface curves step procedures, such 
those Mr. Stevens, would more than minor convenience. 

order illustrate the general applicability these functions, reference 
made graphical method described Woodward, Hon. ASCE, 
and ASCE. Their general equation can written 


which 


The relation between Eqs. and Eq. obvious. 

the graphical solution outlined Messrs. Woodward and Posey, values 
f(y) are obtained for the range depths covered the water surface profile. 
These values f(y) are then plotted ordinates, versus depths abscissas. 
The area under the resulting curve between and the distance between the 
corresponding stations. 

should emphasized that the functions and are general purpose, 
and can applied advantage not only the writer’s solution but the 
Woodward-Posey method and the Husted method well. This point can 
best clarified considering typical problem. 

Example 3.—Water flowing from reservoir into trapezoidal spillway 
channel whose bottom width ft. The slope the spillway floor 
0.0530 and the roughness the lining given 0.025. The rate flow 
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The depth the control section given 4.43 ft. The constant 


(1,000)? 
31,100. 


Values 


100 300 500 1000,10 100 300 500 1000 


The values the functions and are obtained from Fig. which was 
derived from the writer’s the depth factors for side slopes 1.5. 
The surface slopes other depths are obtained similarly. summary 
the calculations required for direct plotting presented Table 13(a). The 
curve the water surface the spillway drawn Fig. 14. 


TABLE FOR EXAMPLE 


Average 
u 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 
4.43 32.1 4.00 29.6 38.9 2.0 
4.00 35.0 45.6 —0.1032 9.68 3.50 45.9 58.6 10.3 
3.50 56.7 71.6 —0.0301 33.2 3.00 77.5 95.8 31.8 
3.00 119.9 —0.0092 108.6 2.75 115.9 139.9 41.5 
2.75 133.6 159.9 —0.0038 261 
2.50 186.6 218.5 0.0000 Summation, 85.6 


order solve this the graphical device suggested 
Messrs. Woodward and Posey, the values f(y) are plotted against the depths 
Fig. 15. The area (A) under the curve equal the distance from the 
control section the station where the depth 4.00 ft; the area (B) equals 

Tables for the values the factors and are available for the cases vertical side walls and side 


and depths from 0.5-ft intervals. The cost these tables $1.00. 


a 
ay 
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the distance from where 4.00 where 3.50 ft, The water sur- 
face profile determined this method superimposed Fig. for com- 
parison. 

For finding the water surface use the Husted method, Eqs. and 
can again advantageous. customary the step methods, the 
values the variables will averaged for each length channel. The 
factors and Eq. therefore, will taken the averages the func- 
tions the two stations that define the increment. The results summarized 
Table 13(b) were obtained applying Eq. 51a finite terms. 


200 280 


2.0 2.75 3.50 4.0 4.50 


Husted and Woodward Posey 


Values 
Values f(y) 


EXAMPLE 


The lengths shown Table are seen practically the same 
those Fig. 15, the latter been computed suggested Messrs. 
Woodward and Posey. 

This typical problem has thus been solved the Husted method, the 
Woodward-Posey graphical procedure, and direct plotting the water 
surface slopes. each these three approaches, the general depth functions 
and have been used reduce the labor computation. The close con- 
formity the resulting water profiles, illustrated Fig. 14, that any 
these solutions provides acceptable accuracy. Whether one another 
method will best solve given case will depend upon the length the profile 
and upon the facilities available. 

These closing statements serve show that Mr. Stevens and the writer are 
not widely separated their lines attack. The two methods stem from the 


200 
160 
| 
\IS | 
4.0 
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same basic relationships and produce practically identical results. each 
case the calculations are started from known depth, and the depth functions 
for each increment are substituted into formula such Eq. The effects 
velocity distribution and eddy loss, predictable, can taken into account 
with equal ease either method. 

Mr. Stevens’ distaste for lengthy formulas and laborious calculations 
shared -by the writer and the multitude. The merits the suggestions 
offered this paper and its discussions will found where there saving 
the time spent the designer. 

Messrs. Bishay and Babbitt have contributed interesting comparison 
several methods use single example involving relatively long water 
surface curve. They have shown the writer’s procedure require probably 
the least volume numerical calculation; but, the same time, they indicate 
appreciable deviation the computed length along the channel for given 
depth. This type comparison misleading because considers only 
relative horizontal distances instead relative depths. this method 
would possible show that two curves approaching near the normal depth 
could differ length percentage approaching infinity, even though the 
actual depths given station were nearly identical. This ludicrous 
example, course, but should illustrate the uncertainties such compari- 
son. 

The practical engineer wants know how deep the stream point 
11,000 down the channel, not how far where the depth 7.0 ft. 
wants know how high build his sidewalls and how much freeboard will 
need. The percentage deviation depths computed here would well 
within the limits usually required the designer. 

Mr. Nelidov has ably stated, the only correct solution integra- 
tion. However, common methods integrating the water surface equation 
involve their own approximations and give accuracy which not always ade- 
quate compensation for the additional labor required. most practical 
cases, ease calculation more attractive than laborious and unnecessary 
refinement. 

the introduction the method direct plotting the water surface 
curve, the tangents were arbitrarily made intersect the midpoints the 
depth increments order attain uniform results for illustration. was 
stated that this procedure can modified practice give more accuracy 
that part the curve where the slope varying problems solved 
the writer through use the unmodified procedure, the discrepancy 
depth has not been greater than 4%, the average being about 2%. more 
accuracy required, smaller depth increments can taken. 

Mr. Pietrkowski states that direct plotting restricted special shape 
cross section. The method fact free any such restriction. Although 
the final equations presented here limited trapezoidal channels, the po- 
tential applications the basic idea are manifold. 

The thought-provoking international responses this paper have en- 
hanced its value 
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